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Abstract 

Evolution of purity in case of factorizable interaction between an open system 
and an environment is investigated. We derive a simple expression for purity de- 
crease at the first instants of evolution (when purity is close to unity). We find that 
purity at very small times is maximal when an initial state of an open system coin- 
cides with one of the eigenstates of the interaction operator, no matter how weak the 
interaction is. On the other hand it is widely known that in general eigenstates of 
the interaction are not pointer states. Therefore the procedure of selecting pointer 
states by purity maximization (known as "purity sieve") should not rely on short- 
time purity behavior. We propose a modification of the purity sieve criterion which 
approximately takes into account purity evolution at longer times. As an example 
of its applicability we recover known results for pointer states of a particle under- 
going quantum Brownian motion; we point out that the criterion is not applicable 
for some other models, however. It is argued that the proposed modified purity 
sieve may be used for selecting pointer states of a particle undergoing decoherence 
through collisions. 

1 Introduction 

Purity sieve p][2] is one of the criteria for selecting pointer states of an open quantum 
system, i.e. the states most robust against decoherence. In this work purity sieve is inves- 
tigated in the case of factorizable interaction between an open system and an environment. 
This case is relevant for the majority of models usually considered in the context of deco- 
herence, for example for quantum Brownian motion model, spin-boson model and some 
spin-spin models. The paper is organized as follows. In the next section we briefly review 
the conventional setting and the basic concepts of the decoherence program, outline the 
role of purity and remind the definition of the purity sieve. In section 2 we analyze the 
short-time behavior of purity. Also a modified (with respect to the canonical purity sieve) 
criterion for selecting pointer states is proposed in this section. In section 3 we introduce 
a model which has the following features: (1) it may be reduced to the "almost factoriz- 
able" form; (2) under certain assumptions decoherence in this model may be regarded as 
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decoherence through collisions; (3) the particular case of this model is quantum Brownian 
motion. We use the proposed modified criterion to find pointer states for this model. 



2 Some basics of decoherence. 

In this section we give a very brief and incomplete sketch of the basic setting of the deco- 
herence program. One may find a number of extensive reviews (e.g. |3j) and pedagogical 
introductions to the topic (e.g. [1]). In particular, a profound list of references to the 
original works is given in [3]. 

Consider a System (with Hilbert space S), which interacts with an Environment (with 
Hilbert space 8). The total Hamiltonian of the combined closed system (with Hilbert 
space Ti. = S <8> 8) reads 

H = H s + H £ +V, (1) 

H s , H £ and V being self-Hamiltonians and interaction Hamiltonian correspondingly. The 
usage of superscripts S and 8 here and in what follows is self-explanatory. 

The state of the combined system 7i may be described by a state vector \l/ or by a 
density matrix p. The states of the System S and the Environment 8 are described by 
the reduced density matrices 

p s = tr £ p, p £ = tr s p. (2) 
Density matrices evolve according to the master equations: 

0) 

±pS = -i[H s ,p s ]-it u [V,p}. (4) 

According to the Schmidt theorem, at any time p s and p £ may be diagonalized in 
some basis: 

P S (t) = PMPfit), (5) 

i 

P £ (t) = E Mt)Pf(t)- (6) 

i 

Here {Pf{t)} and {Pf(t)} are systems of mutually orthogonal projectors, and {pi{t)} is 
a set of non- negative numbers (probabilities). 

An essential ingredient of the decoherence program is the following 
Decoherence Hypothesis: for generic System, sufficiently large Environment, generic 
interaction and (almost) any initial conditions, the projectors {Pf(t)} almost always be- 
long to some set of fixed projectors, which correspond to the so-called pointer states. 

In other words, the reduced density matrix is diagonal almost always in the basis of 
pointer states. Exceptional periods of time when the diagonalization basis substantially 
deviates from the pointer basis are short and rare. If the initial state of the System is 
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far from pointer, the first moments of the evolution form such an exception: the System 
decohere in the mixture of the pointer states during this period. 

One of the main problems of the decoherence theory is to find pointer states, given 
a Hamiltonian H of the combined system 7i. One of the helpful tools for this is purity 
sieve PQ,[2]- In order to introduce purity let us first note that there is a simple method to 
obtain the probability p max , which is by definition the largest one among all Pi, and the 
corresponding projector P^ ax , given the reduced density matrix p s : 

PL* = I** [iP S ) k /^s(P S ) k ] , (7) 

k— >oo 

JW = tr s Ptj S = lim [tT S (p s ) k+1 /tr s (p s ) k ] . (8) 
In practice, approximate expressions for -P^ax an d Pmax may be used: 

PL, * (p S ) k /tr S (p S ) k (9) 
2W * tT S (p S ) k+1 /trs(p S ) k . (10) 

Here the larger is k, the better is approximation. 

Assume now that at t — the state of the System is pure, p 5 (0) =| ijj s )(ip s \ . Then 
p max initially equals 1, and stays close to 1 during some time. In this case even k = 1 
appears to be a good approximation, and one gets 

Pmax ^ tl S (p S ) 2 . (11) 

The quantity in the r.h.s. of the above equation is called purity. We denote it by : 

V = tr 5 (p 5 ) 2 - (12) 

Sometimes —^3 is called " linear entropy" . 

For the pure state p s =\ tp s )(ip s \ purity equals 1, and it decreases in the course 
of evolution. The rate of the decrease differs for different initial states. It is natural 
to assume (in the spirit of the Decoherence Hypothesis) that for pointer states this rate 
should be small. The small rate of purity loss implies that the corresponding state is 
stable against decoherence. This reasoning suggests that in order to find pointer states 
one should minimize the purity decrease rate with respect to the initial pure states. This 
procedure is called purity sieve. 

To summarize, purity sieve is a procedure for selecting pointer states of the System. 
It comes to the following two steps: 

1. calculating purity ^3 as a functional of the initial state of the System ip s , a functional 
of the initial state of the Environment p £ , and a function of time £; 

2. maximizing with respect to ip s for some fixed (reasonable) t and p . 

The states which provide the local maxima should be the pointer states. The procedure 
makes sense if the result weakly depends on p and t, provided t is greater than some 
characteristic time, called decoherence time. 

lr This implies that the combined closed system starts its evolution from the initial product state 
p = \tl) s ){^ s \ ®p £ . 
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3 Purity sieve for factorizable interactions 



A number of models, investigated in the context of decoherence studies, have factorizable 
interaction Hamiltonian (see e.g. [2], [5j-[8j): 

V = SE. (13) 

Such form of interaction allows to simplify the master equation (jl]): 

j t p s = -i[H s ,p s }-i[S,tr £ (Ep)}. (14) 

Using this expression, we calculate the first two derivatives of purity at t = : 

# = 0, (15) 

% = -45E 2 5S 2 . (16) 

Here the dispersions squared are calculated for the initial states of the System and the 
Environment: 

5E 2 = tr £ E 2 p £ - (ti £ Ep £ ) 2 , (17) 
5S 2 = | S 2 | - I S | i) s ) 2 . (18) 

The first derivative of purity vanishes at t — merely because the evolution starts 
from the state with maximal purity, *p = 1. 

The formula for the second derivative, ffTB"]) . is a result of the straightforward calcula- 
tion; it appears to be less lengthy if one uses the interaction representation instead of the 
Schroedinger representation. 

The short-time behavior of purity is quadratical in time and reads 

*P ~ 1 - 25E 2 5S 2 t 2 . (19) 

A noticeable fact is that in the case of factorizable interactions the second derivative at 
t — also factorizes. As a consequence, as long as eq. (1191) holds, purity decrease is minimal 
for eigenstates of the interaction operator S, which give SS 2 = and ^3=1 — 0(t 3 ). Does 
this mean that purity sieve selects the eigenstates of S to be the pointer states of the 
System? In general, the answer is negative. For example, if the interaction is weak 
compared to the self-Hamiltonian of the System, pointer states are known to be the 
eigenstates of the self-Hamiltonian H s [Sj. 

To illustrate this situation let us consider a simple central-spin model as an example: 

N N 

# = !**+£^ +aT "I>£- (20) 
i=i i=i 

Here a x and a z correspond to a central spin, a l x and a\ correspond to iV bath spins, con- 
stants to and e characterize the strengths of self-Hamiltonian and interaction correspond- 
ingly. ^3(t) is calculated numerically for two initial directions of the central spin, along Oz 
axis, \0z) (the eigenstate of the self-Hamiltonian H s = %cr z ), and along Cte axis, \0x) (the 
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purity 




Figure 1: Purity as a function of time for the central spin model ( 1201) and for two different 
initial states, \0z) (thick green line) and \0x) (thin red line). The parameters of the model 
are N = 6, ui — 1, e = 0.1. The initial state of the environment is a mixture of all possible 
states with equal weights (p £ = 1/2 N ). 

eigenstate of the interaction operator S = a x ). The result for N = 6, u> = 1, e = 0.1 is 
shown in Fig.l. One can see that at first moments the eigenstate of the interaction opera- 
tor, | Ox), indeed "wins" the competition, providing minimal purity decrease. However, at 
longer times the situation is not so clear. Purity evolution is no longer described by the 
quadratic law ffT9l) ; purity sieve does not unambiguously select any state. We conclude 
that eq. (1191) correctly describes short-time behavior of the considered model, but purity 
sieve is hardly applicable in this case. 

Another example is a particle which undergoes decoherence through collisions with 
other particles. The interaction hamiltonian in this case is a function of distance between 
particles. According to eq. (fT9l) sharply localized states (coordinate delta- functions) max- 
imize purity at small times. However, such wave functions are known to instantly spread 
over large distances. This implies an abrupt decrease of purity at slightly longer times. 
Therefore, in this case one also can not use eq. (|T9l) to select pointer states. 

The examples considered show that maximizing short-time expression ( JT9l) does not, 
in general, provide correct pointer states. However, the later example gives a hint how 
one can try to upgrade a selection criterion. The key quantity in eq. (1191) is the dispersion 
squared of the interaction operator, 5S 2 . In order to maximize purity we minimized SS 2 , 
calculated for the initial state i/j s . A straightforward generalization is to try to calculate 
the dispersion of the interaction operator for longer times, to average it and then - to 
minimize it. Thus, a modified criterion for selecting pointer states may be proposed, 
which comes to the minimization of 

f dt'5S 2 (t'), (21) 
Jo 

where 

5S 2 = (^ s (t) | S 2 | V 5 (*)> - (^ S (t) | S | ^(t)) 2 , (22) 

2 The factor SE 2 is expected to take some typical value for large environments; it should not affect the 
result. 
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and ip s (t) is an evolving state of the System (which is assumed to be approximately pure). 
It satisfies a Shroedinger-like equation 

ij t ^ S (t) = (H s + nS)ip s (t), k = tT £ Ep £ . (23) 

This criterion appears to work satisfactory in the case when the interaction operator 
is a coordinate operator, S — x. This may be verified by considering specific examples 
(see the next section). On the other hand, one can easily check that it does not work, 
for example, in spin-spin model (T2T)]) . The reason for such diverse behavior seems to be 
connected with the singularity or non-singularity of the interaction operator. However, 
this question is far from being clearly understood. It will be discussed elsewhere. 



4 Purity sieve for quantum Brownian motion and for 
decoherence through collisions 

In this section we give an example of implementation of the the proposed modified selection 
criterion. 

The following Hamiltonian is considered: 



2 i N 2 i N N 

6= 2M + 2 MfiV + + 2 mUj2xl) + E U ^ Xk ~ x ) + E U ^ Xk ~ ^ ( 24 ) 

k=l k=l k>j 



It describes one particle of mass M and iV particles with masses equal to m in the 
oscillator potential wells. All particles interact with each other through potentials XJ\ and 
IJ%. If this potentials are collision potentials, the model describes decoherence through 
collisions. If the potentials are quadratic, this is a quantum Brownian motion (QBM) 
model. 

Let us rewrite this Hamiltonian in an "almost factorizable" form. For this we use the 
following transformation of variables ([21 

x = x 

•&k %k X 

P = p + ^2Pk 
Pk = Pk- 

In new variables the Hamiltonian reads 



H = H s + H £ + V, 

H s = -C- + -(Mfi 2 + Nmuj 2 )5: 2 (25) 

2M 2 V ; y ' 

N _ 2 _ \2 i N N 

k& = E(£ + -^Sr + \ muj2il) + E u &») + E - ^ ( 26 ) 

k=l ' ' k=l k>j 

V = — hmw x2_^x k . (27) 
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One can see that the interaction Hamiltonian is a sum of two factorizable terms now. 
We could make some further assumptions in order to make one of the terms in (j27p 
negligible and get rid of it. However, it appears to be unnecessary. Indeed, let us calculate 
the integral (1211) . which is subject to minimization, separately for both terms: 

F dt'8f{t') = l/2(M 2 Q 2 5x(0) 2 + Sp(0) 2 )T, 
Jo 

F dt'5x 2 (t r ) = l/2(5x(0) 2 + 5p(0) 2 /M 2 n 2 )T, 
Jo 

where 

fi 2 = & + (28) 

and integration is performed over the oscillator period T. We see that due to p — x 
symmetry in the oscillator Hamiltonian two contributions appear to coincide up to a 
constant factor. One can prove that the interference term vanishes. Minimizing the 
obtained expression one gets the pointer states to be those with 

5x 5p = 1/2, (29) 

8x 2 = (2Mny\ (30) 

Thus pointer states are minimum uncertainty states with the coordinate uncertainty 
defined by (1301) . 

Previously this result was obtained for the QBM [2J. From the mathematical point 
of view calculations in [2] coincide with those in the present paper, although a different 
approach, based on the approximate master equation, was used in [2]. 

Approaches, which are also very different from the presented one, are usually used to 
study decoherence through collisions (see, e.g. [Tj , [TO] ) - These studies normally predict 
the pointer states to be phase-space localized states (as demanded by the quantum-to 
classical correspondence principle), which agrees with our result. 

Although our result for the model (l2~4"l) seems to be in general reasonable, one should 
be cautious when applying it to specific forms of Hamiltonian with certain functions 
U\ and U2- For example, if the interaction vanishes, U\ =0, the very conception of 
pointer states makes no sense, but our selection criterion does not "feel" it and yet gives 
minimum- uncertainty states as pointer states. As it was mentioned above, the applica- 
bility conditions of the proposed criterion are not entirely explored. 
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